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ABSTRACT 

We perform a combined likelihood analysis of the latest cosmic microwave background 
anisotropy data and distant Type la Supernova data of Perlmutter et al. (1998a). Our 
analysis is restricted to cosmological models where structure forms from adiabatic ini- 
tial fluctuations characterised by a power-law spectrum with negligible tensor compo- 
nent. Marginalizing over other parameters, our best fit solution g 

and n A = 0.63^23 (95% confidence errors) for the cosmic densities contributed by 
matter and a cosmological constant respectively. The results therefore strongly favour 
a nearly spatially flat Universe with a non-zero cosmological constant. 



1 INTRODUCTION 



The aim of this paper is to constrain the geometry of the 
Universe by combining results from the cosmic microwave 
background anisotropies (CMB) with those from distant 
Type la supernovae. Recently, there have been a number 
of analyses of CMB anisotropies aimed at constraining cos- 
mological parameters (Hancock et al. 1998, Bond and JafFe 
1997, Lineweaver and Barbosa 1998a,b, Webster et al. 1998). 
These papers show that the CMB already provides useful 
constraints on adiabatic inflationary models. The results de- 
scribed in this paper extend these analyses to a wider pa- 
rameter set including closed Universes. 

However, observations of CMB anisotropies alone can- 
not determine the geometry of the Universe unambiguously 
(Bond, Efstathiou and Tegmark 1997, Zaldariaga, Spergel 
and Seljak 1997). This is because two models with identi- 
cal fluctuation spectra and matter content will have nearly 
identical CMB power spectra if they have the same angular 
diameter distance to the last scattering surface. This ge- 
ometrical degeneracy can be broken for extreme values of 
the cosmological parameters £Ia and Q m by the inhomoge- 
neous Sachs-Wolfe effect at low multipoles (see Efstathiou 
and Bond 1998), but for plausible parameters the degener- 
acy is nearly exact. ^ 

A number of authors (White 1998, Tegmark et al. 1998, 
Efstathiou and Bond 1998) have shown that the magnitude- 
redshift relation for distant Type la supernovae (SN) pro- 
vides nearly orthogonal constraints in the QA-Qm plane to 
those derived from the CMB. The combination of CMB and 
SN data can thus provide tight constraints on the geome- 



* Gravitational lensing of the CMB can break this degeneracy 
however (see Seljak and Zaldarriaga 1998). 



try of the Universe. This has been demonstrated by White 
(1998), Garnavich et al. (1998) and Lineweaver (1998) us- 
ing the SN data of Perlmutter et al. (1997, 1998b) and Riess 
et al. (1998). In this paper, we use the larger sample of 42 
high-redshift SN of Perlmutter et al. (1998a, hereafter P98, 
1998b). The likelihood analysis of the CMB observations is 
described in Section 2. Section 3 describes the likelihood 
analysis of the SN data and the results from the combined 
data set are described in Sections 4 and 5. 



2 ANISOTROPIES OF THE COSMIC 
MICROWAVE BACKGROUND 

The analysis presented here is similar to that described by 
Hancock et al. (1998) and we refer the reader to this pa- 
per for technical details. Similar analyses and compilations 
of observations are discussed by Lineweaver (1998, and ref- 
erences therein) and by Bond, Jaffe and Knox (1998). A 
window function Wi for each experiment is used to convert 
the observed level of anisotropy to flat bandpower estimates 
(ATi/T)±a centred on the effective multipole Z e g (defined as 
the half-power point of the window function). The resulting 
CMB data points are plotted in Figure 1, together with their 
68 per cent confidence limits. These confidence limits have 
been obtained using likelihood analyses and hence incorpo- 
rate uncertainties due to random errors, sampling variance 
and cosmic variance. 

The data points in Figure 1 are identical to those given 
in Webster et al. (1998) except that we have added the three 
QMAP points (Ka and Q bands for flights 1 and 2, Devlin et 
al. 1998, Herbig et al. 1998, De Oliveira-Costa et al. 1998). 

We fit the CMB data points to adiabatic CDM models 
specified by the following parameters: (1) amplitude Qio of 
the I = 10 multipole defined as in Lineweaver (1998); (2) the 
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Figure 1. CMB bandpower anisotropy estimates for various ex- 
periments, as described in the text. The solid line shows the 
best fitting adiabatic CDM model with parameters Qio = L14, 
n = 1.08, lu c = 0.36, u>b = 0.03 and Doppler peak location pa- 
rameter 7o = 1.18. The dotted lines show the best fitting curves 
with location parameter fixed at 70 = 0.8 and 70 = 1.5 (approx- 
imately the 2<r allowed range for 7/3 . 



density parameters fij. and (3) J1a; (4) the scalar spectral 
index n s ; (5) the physical density in cold dark matter ui c = 
Q, c h 2 [j (6) the physical density in baryons Ub = 0.bh 2 . We 
ignore tensor modes in this analysis.^] 

The motivation for choosing this set of variables is ex- 
plained in Efstathiou and Bond (1998, hereafter EB98). The 
physical densities uj c and uib and the radiation density de- 
termine the sound speed at the time of recombination. The 
geometry of the universe is specified by the parameters £lk 
and S1a, and the Hubble constant enters as an auxiliary pa- 
rameter 



h 



(1 - O fe - fi A ) 



1/2 



(1) 



We have fitted the observations to the theoretical mod- 
els using two methods. Firstly, we compute a grid of theoreti- 
cal power spectra in the parameters Qk, ^a, n s ,u) c , using the 
CMBFAST code of Seljak and Zaldarriaga (1996) with uj b 
constrained to 0.019, the value inferred from primordial nu- 
cleosynthesis and the deuterium abundances measured from 
quasar spectra (u>b = 0.019 ± 0.001, see Buries and Tytler 
1998a,b). In its present form, the CMBFAST code is re- 
stricted to open and spatially flat models (fife > 0) and so 
we have adopted a second, approximate, technique to ex- 
tend the theoretical predictions to closed models. This is 
based on a semi-analytic fitting formula for the CMB power 
spectrum of Qk — models, which is a generalization of the 
fitting formula of equation (24) in EB98 and provides accu- 



t Here h is the Hubble constant in units of lOOkms 1 Mpc 1 . 
T For reasonable amplitudes of a tensor mode, the effect on the 
position of the first Doppler peak is small. A small tensor mode 
will therefore have little effect on the cosmological parameters £2 a 
and Q m since these are determined primarily by the position of 
the first Doppler peak. 



Figure 2. Likelihood contours in the fl\ — O m plane derived 
from the CMB data points shown in Figure 1. The contours are 
plotted where —21nC/C m ax is equal to 2.29, 6.16 and 11.83, cor- 
responding approximately to 1, 2 and 3<r confidence contours for a 
Gaussian likelihood function. The solid lines show the marginal- 
ized likelihood contours derived from the CMBFAST computa- 
tions for Ub = 0.019. The dotted lines extending into the < 
region show the equivalent contours derived from the fitting func- 
tion approach described in the text. For the dotted contours, we 
have marginalized over Ub, although this has very little effect. 

rate fits to the first three Doppler peaks of the CMB power 
spectrum. The fitting formula includes the dependences on 
n s , uJb, ijj c and Qa and is typically accurate to better than 5 
percent. For models with non-zero values of Qk, we use the 
scaling relation C(t!) — + C(£jd), where 713 is a 'location' 
parameter, 



7d 



e D (n k = o,si A = o)' 



(2) 



where ^(fij^O/v) is the location of the first Doppler peak 
given by equation (22) of EB98 generalized to models with 
Qk < 0. The location parameter thus measures the positions 
of the Doppler peaks relative to those of a spatially flat 
model with zero cosmological constant. The approximate 
formula does not include the inhomogeneous Sachs- Wolfe ef- 
fect (see e.g. Bond 1996) which affects low multipoles if £lk 
and Q,a are non-zero. However, the inhomogeneous Sachs- 
Wolfe effect is a weak discriminator of cosmological mod- 
els (see EB98) and, with the data shown in Figure 1, the 
constraints on cosmological models are set primarily by the 
location of the first Doppler peak. A similar approximate 
technique, using rescaling of power spectra computed with 
CMBFAST, is described by Tegmark (1998). 

Figure 1 shows the best fitting CMB power spectrum, 
together with the best fitting curves with 713 = 0.8 and 
1.5 (spanning the 2a allowed range of after marginaliza- 
tion over other parameters). The present CMB data points 
evidently favour models with Qk ~ 0. The large number 
of data points at £ ~ 100 set quite strong constraints on 
closed models, but the lack of data points at £ £ 300 leads 
to weaker constraints on open models. The CMB data do 
not yet allow strong constraints on the parameters io b and 
uj c , hence the constraints on Hubble constant (equation 1) 
are also extremely weak. 

Figure 2 shows the CMB constraints in the Q m — Q,a 
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plane. Here, we have marginalized over n s , Qio, u c (and 
uib for the likelihoods computed from the fitting formula) 
assuming uniform prior distributions in these parameters. 
The marginalized likelihood depends slightly on the range 
adopted for u> c ; in Figure 2 we assume a uniform prior dis- 
tribution over the range 0.05 < u> c < 0.5. The marginalized 
likelihood function is insensitive to the ranges and priors 
adopted for the other parameters and is insensitive to LUb- 

The results of this analysis are similar to those of 
Lineweaver (1998) and Tegmark (1998) but differ in detail. 
The main difference is in the way that we marginalize over 
the likelihood function. We have assumed a uniform prior 
distribution in each parameter and performed direct inte- 
grations over the full likelihood function. Lineweaver and 
Tegmark 'marginalize' over a subset of parameters by fix- 
ing them to their maximum likelihood values. Our approach 
provides more robust errors, though the marginalized likeli- 
hood function for poorly determined parameters will depend 
on the choice of prior (usually weakly). In the limit that 
the likelihood function is Gaussian, the two approaches are 
equivalent. However, for non-Gaussian likelihoods (as is the 
case with the present CMB data) the approach adopted by 
Lineweaver and Tegmark can give misleadingly small errors 
on some parameters. 



3 MAGNITUDE- REDSHIFT RELATION FOR 
DISTANT SUPERNOVAE 

We use the sample of 42 high redshift (0.18 < z < 0.83) 
supernovae of P98, supplemented with 18 low redshift (z < 
0.1) Type la supernovae from the Calan/Tololo Supernova 
Survey (Hamuy et al. 1996). For each supernova, P98 com- 
puted a peak magnitude in the B band tub, corrected 
for Galactic extinction and a 'stretch parameter' s that 
stretches the time axis of a template Type la lightcurve to 
match the observed lightcurve (see Perlmutter et al. 1995, 
1997). 

P98 provide a comprehensive analysis of the constraints 
on Q m and J1a derived from the SN magnitude-redshift re- 
lation and of the effects of excluding various outlying SN, 
including or excluding corrections for the lightcurve width- 
luminosity relation, host galaxy extinction, etc. P98 show 
that the likelihood function in the fi m -fiA plane is remark- 
ably stable to such changes. We do not repeat this analysis 
here, but instead concentrate on the analysis of the super- 
novae used in the 'primary fit' of P98 (their fit C) which ex- 
cludes 4 high redshift objects. These are, SN 19970, 1996cg 
and 1996cn, which are very likely reddened by their host 
galaxies and so are fainter than the best fitting magnitude 
redshift relation, and SN 1994H which is not spectroscopi- 
cally confirmed as a Type la SN and lies brighter than the 
best-fit relation. In agreement with the results of P98, none 
of our conclusions change significantly if we include these 
supernovae. 

We define a corrected peak magnitude m c g TT for the 
lightcurve width-luminosity effect 

ma" = m B + a(s - 1), (3) 

where s is the measured stretch factor and a is a constant 
to be determined. These corrected magnitudes are compared 
to the predicted magnitudes 
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Table 1: Maximum Likelihood Parameters 



M B 


a 




n A 


sample 


-3.45 


1.33 


0.54 


1.09 


combined 


-3.69 


1.41 


0.75 


1.86 


high z 


-3.37 


1.32 






low z 


= Mb 4 


- 51o g r>i 




,Oa), 





where Mb is related to the corrected absoluted magnitude 
Mb by Mb = Mb - 51og# + 25, and V L = d L + 51og// 
is the Hubble-constant-free luminosity distance defined by 
P98. To compute the luminosity distance, we ignore gravita- 
tional lensing and use the standard expression for a Universe 
with uniform density (see e.g. Peebles 1993), 

d L (z,n m ,n A ) = |^ sin fc [\Qk\ 1/2 x(z,fl m ,n A )] , 
/* z dz 

x(z,n m ,n A ) = [Qm{1 + z , )3 + Qkil + zl)2 + a A] i/2 ^ 

where Qt = 1 — f!m — and sin^ = sinh if fife > and 
sink = sin for flj. < 0. We assume a constant cosmological 
constant here rather than an arbitrary equation of state as 
might arise with scalar fields that become important at late 
times (Ratra and Peebles 1988, Caldwell, Dave and Stein- 
hardt 1998). The biases in the magnitude-redshift relation 
arising from gravitational lensing should be negligible for 
CDM-like models unless a large fraction of the dark matter 
is in compact objects (e.g. Wambsganss, Cen and Ostriker 
1998). Even in the latter case, P98 show that the biases are 
relatively small for the low matter densities favoured by the 
SN data. 

With the above assumptions, we perform a four parame- 
ter (fl m , Q A , Mb and a) likelihood analysis assuming Gaus- 
sian errors on m™ r consisting of three terms, 

(Am C B rr ) 2 = Am| + a 2 As 2 + Amf ntrinsic , (6) 

where Ams and As are the measurement errors in ms and 
s and Amintrinsic is the instrinsic dispersion in tub deter- 
mined to be 0.18 mag from the Calan/Tololo sample. (The 
maximimum likelihood parameters are extremely insensitive 
to Amintrinsic-) This analysis differs from that in P98 in that 
we include the dependence of the magnitude errors on the 
parameter a self-consistently in the likelihood analysis via 
equation (6) . As we will see below this has little effect on the 
likelihood constraints on fl A and Q, m (Figure 3), which are 
in good agreement with the results presented in P98. How- 
ever, our analysis allows an important test of the intrinsic 
properties of high and low redshift supernovae (see Figure 
4). 

Excluding the four SN as described above, we analyse 
the combined high redshift P98 + Calan/Tololo samples (de- 
noted 'combined'), and the two samples separately (denoted 
'high z' and 'low z' respectively). The magnitude-redshift 
relation of the low z sample is relatively insensitive to cos- 
mology, hence we fix Q m = 1 and £l A — in the likelihood 
analysis. Table 1 lists the parameters that maximize the 
likelihood for these samples. 

Figure 3 shows the likelihood function for the combined 
sample in the Q m -Q, A plane after marginalization over the 
parameters M b and a assuming uniform prior distributions 
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Figure 3. Likelihood contours (1, 2 and 3cr, as denned in the 
caption to Figure 2) in the Q^-Qm plane derived from the SN 
magnitudc-rcdshift relation after marginalization over the param- 
eters Mb and a. The cross shows the maximum likelihood so- 
lution. Bouncing universes have parameters within the hatched 
region. 



in these variables. We use this likelihood function in Sec- 
tion 3 when we combine the SN sample with parameters de- 
rived from the CMB. As P98 show (and we have confirmed) 
changes in the analysis, e.g. omitting outliers, stretch cor- 
rection, reddened objects, usually shifts the error ellipses by 
much less than the width of the la contour. 

An important consistency check is the agreement be- 
tween the intrinsic properties of the high and low redshift 
SN. This is illustrated in Figure 4 which shows likelihood 
contours for the parameters Mb and a for the low z and 
high z samples (marginalized over Qa and Q m for the high z 
sample) . This diagram shows that the low z and high z sam- 
ples have the same lightcurve width-luminosity relation and 
consistent peak absoluted magnitudes M b ■ (Note, however, 
that the contours for the high z sample become highly elon- 
gated in the M b direction because this parameter correlates 
strongly with Qa and Q m ). We have repeated the likelihood 
analysis including the intrinsic magnitude scatter Ami ntr insic 
as a free parameter. For the low-z sample, the likelihood 
gives a broad distribution peaked at ZXTTlintrinsic = 0.18, 
consistent with the rms residual of 0.16 magnitude of the 
Calan/Tololo points around the best fit solution. The com- 
bined sample gives Ami Iltr insic = 0.20, slightly higher but 
consistent with the low-z sample. (In fact, Amintrinsic drops 
to 0.18 mag if we remove three outliers, SN92bi, SN95as and 
SN97K). 

For low redshift supernovae, the lightcurve widths and 
luminosities are known to correlate with the morphological 
type of the host galaxies (Hamuy et al. 1996). Thus the lack 
of any significant difference between the high and low red- 
shift samples in Figure 2 provides a constraint on evolution- 
ary corrections associated with systematic changes in the 
host galaxies with redshift (e.g. metallicity, morphological 
mix). Figure 4 therefore indicates no detectable evolution- 
ary trends in the underlying supernova population used in 
this analysis. 

Table 1 shows that the high z sample alone yields val- 




Figure 4. Likelihood contours (1, 2 and 3cr) in the Mb~ & plane. 
The solid lines show contours for the Calan/Tololo sample and 
the dotted lines show contours for the high redshift P98 sample 
after marginalization over the cosmological parameters Qa an d 
f2 m . 



ues for Q m and Qa that are similar to those of the combined 
sample. However, the error contours for the high z sample 
are much larger than those shown in Figure 1 for the com- 
bined sample so that an Q m — 1, Qa = universe lies within 
the 2a contour. The narrowness of the likelihood contours 
shown in Figure 3 therefore rely on combining the P98 data 
with the Calan/Tololo sample, though the general trends are 
evident in the high-z sample alone. 

In summary, the P98 and Calan/Tololo SN strongly 
favour a universe with 0.78fi m — 0.62^a ~ —0.25 ± 0.13. 
This is consistent with the analysis of P98 (who find 0.8O m — 
0.6S7a ~ —0.2 ± 0.1) and with the analysis of a sample of 
16 SN at z > 0.16 (14 of which are independent of the P98 
high redshift sample) combined with 34 low z SN from the 
Calan/Tololo and CfA samples (Riess et al. 1998, Garnavich 
et al. 1998). Furthermore, the lightcurve width-luminosity 
relation of the high z SN is consistent with that of the nearby 
sample suggesting that these objects have similar intrinsic 
properties. If an evolutionary effect is causing a systematic 
error in the magnitude-redshift relation of the distant sam- 
ple, then it must be so as to preserve the lightcurve width- 
luminosity relation. 



4 COMBINING THE SUPERNOVAE AND 
CMB LIKELIHOODS 

4.1 Combined likelihoods 

The combined likelihood obtained by multiplying the SN 
and CMB likelihoods are shown in Figure 5. As in Figure 
2, the solid lines show the likelihood computed with CMB- 
FAST and the dashed lines show the approximate technique 
indicating how the contours extend into the Qk < regime. 
The likelihood peaks at Q m = 0.25 and Qa = 0.63. The com- 
bined likelihood thus strongly favours a nearly spatially flat 
Universe with a low matter density and high cosmological 
constant. In fact, the 2a ellipse in Figure 5 extends over the 
range Q m « 0.12, Qa = 0.84, to Q m « 0.49 and Qa = 0.51. 
A high value of Qa is suggested by the SN data alone, and is 
required if we impose the constraint Q k = (see Figure 7 of 
P98; Figure 6 Riess et al. 1998). However, from the SN data 
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Figure 5. Likelihood contours derived by combining the super- 
novae likelihood function shown in Figure 3 with the CMB like- 
lihood function shown in Figure 2. As in Figure 2, solid contours 
show 1, 2 and 3<x confidence intervals computed using CMBFAST. 
The dashed contours show 2 and 3<r contours computed using the 
approximate CMB fitting technique (for clarity we do not plot the 
la contour). The combined CMB+SN likelihood function peaks 
at Q m = 0.25 and Qa = 0.63 

alone we cannot rule out an open Universe with a low mat- 
ter density Q m £ 0.1 and zero cosmological constant. Since 
the CMB data favour a universe with Qk = 0, the combined 
SN+CMB data require a non-zero cosmological constant at a 
high significance level. This is the main result of this paper. 

Figure 6 provides another illustration of how the combi- 
nation of CMB and SN data dramatically improve the con- 
straints on Qa and Q m . Here we have plotted the likelihood 
functions marginalized over all other parameters except Q m 
(Figure 6a) and fi A (Figure 6b) for the SN and CMB data 
alone and for the combined data sets. For a Gaussian like- 
lihood function, the 95% confidence region is delineated by 
£/£max > 0.146 and so we can see from Figure 6 that the 
constraints on Q m and Qa from the SN and CMB data alone 
are extremely weak. However, for the combined data sets we 
find Q m = 0.25±°;^ and Q A = 0.63tg;^ at the 95% confi- 
dence level. 



4.2 Systematic errors 

Systematic errors in the analysis of the SN data are dis- 
cussed at length by P98 and by Riess et al. 1998. No sys- 
tematic error has yet been identified that could produce a 
large downward shift of the SN likelihood contours in Figure 
3. Possible sources of systematic error include differences in 
the reddening caused by the host galaxies, and evolution- 
ary (e.g. metallicity dependent) corrections to the SN abso- 
lute magnitude-redshift relation. Internal reddening for this 
sample is discussed in detail by P98. They find no excess 
reddening of most of the distant SN when compared to the 
Calan/Tololo sample, although there are a small number of 
possibly reddened SN (three of which have been excluded in 
this analysis). The inclusion or exclusion of these reddened 
objects does not significantly affect Figure 5. Grey extinc- 
tion is much harder to rule out, but may not be physically 
well motivated. 



a n. 

m A 

Figure 6. Marginalized likelihood functions plotted as function 
of fi m (left hand panel) and Qa (right hand panel). The dotted 
lines are for the SN data alone (Figure 3), dashed lines for the 
CMB data alone (Figure 2) and solid lines for the combined data- 
sets (Figure 5). 

Possible evolutionary effects are difficult to check. How- 
ever, the results summarized in Figure 4 show that the high 
and low redshift SN sample have statistically indistinguish- 
able internal properties, i.e. they have consistent lightcurve 
width-luminosity relations and (within rather large errors) 
consistent peak absolute magnitudes Mb- With a larger 
sample of SN it should be possible to refine this test and to 
test for differences in the distribution of lightcurve shapes 
with redshift. Another consistency check would be provided 
by extending the SN sample to z > 1. As Goobar and 
Perlmutter (1995) have discussed, the degeneracy of the 
magnitude-redshift relation in the Q m — Qa plane can be 
broken by a sample of SN spanning a sufficiently wide range 
of redshifts. 

4.3 The best fit Universe 

If systematic errors are indeed small, the combined CMB 
and SN data strongly favour a near-spatially flat Universe 
with Q m m 0.25 and Qa ~ 0.63. These values are close to 
those favoured by a number of other arguments, which we 
summarize briefly below (see also e.g. Ostriker and Stein- 
hardt 1995, P98, and references therein). 
Age & Hubble constant: For our best fit cosmology, the age 
of the Universe is 14.6(/z/0.65) _1 Gyr (in agreement with 
P97). This is compatible with recent estimates of 11.5 ± 1.3 
Gyr for the ages of the oldest globular clusters (see Chaboyer 
1998) and with recent values of Ho derived from Type la su- 
pernovae and Cepheid distances, which fall within the range 
H = 65 ± 10 kms^Mpc -1 (e.g. Freedman et al. 1998). 
Large-scale structure: Observations of large-scale structure 
(see e.g. Efstathiou 1996 for a review) are consistent with 
scale-invariant adiabatic cold dark matter universes if T w 
Q m h « 0.2-0.3. This is broadly consistent with the analysis 
presented here and with the analysis of combined CMB and 
IRAS galaxy data presented by Webster et al. (1998). 
Baryon abundance in clusters: Consistency between primor- 
dial nucleosynthesis (u> b « 0.019) and the ratio of baryons 
in clusters to total cluster mass (/{, f» 0.06/i~ 3 / 2 , sec Evrard 
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1998 and references therein) requires a low matter density, 
fi m ~ 0.26(/i/0.65) _1//2 , consistent with the best fit solution 
of Figure 5. 



5 CONCLUSIONS 

• We have applied an approximate formula for the CMB 
power spectrum that can be used to constrain a wide set 
of cosmological parameters, including closed universes, by 
fitting to the CMB anisotropy data. The results agree well 
with those derived from exact computations using the CMB- 
FAST code. 

• In our analysis we perform a proper marginalization over 
parameters, assuming uniform prior distributions, to derive 
constraints in the f2 m -fiA plane and on fl m and Qa sepa- 
rately. 

• Current CMB anisotropy data provide strong constraints 
on the position of the first Doppler peak and favour a spa- 
tially flat Universe. 

• A likelihood analysis of the SN data provides robust con- 
straints on Q m and JIa consistent with those derived by 
P98. For a spatially flat Universe, the SN data require a 
non-zero cosmological constant at a high level of significance 
(Q A ~ 0.5 at 95 % confidence). 

• The lightcurve width-luminosity relation for high redshift 
and low redshift SN are statistically indistinguishable, con- 
sistent with no evolution of the SN population. 

• The combination of CMB and SN data thus provide strong 
constraints on fl m and Q,a favouring values of tt m ~ 0.25 
and Q,a ~ 0.63. If there are no significant systematic effects 
in the SN data (e.g. grey dust, evolution) then we are forced 
into accepting a cosmological constant or a "quintessence" - 
like component of the Universe (Caldwell et al. 1998, Gar- 
navich et al. 1998). 

• There are a number of independent lines of argument, e.g. 
the age of the Universe, large-scale clustering of galaxies 
and the baryon content of clusters, to support the best-fit 
parameters derived in this paper. 
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